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The unbound nucleus 12 Li is evaluated by using the multistep shell model in the complex energy 
plane assuming that the spectrum is determined by the motion of three neutrons outside the 9 Li 
core. It is found that the ground state of this system consists of an antibound l/2 + state and that 
only this and a l/2~ and a 5/2 + excited states are physically meaningful resonances. 
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I. INTRODUCTION 

The study of halo nuclei is one of the main subjects of 
research in nuclear physics at present. Many theoretical 
predictions on halo, superhalo and antihalo nuclei have 
been advanced in recent years Most of these cal- 

culations correspond to nuclei very far from the stability 
line. They are mainly thought as a guide for experiments 
to be performed in coming facilities. The general feature 
found in these calculations is that a necessary condition 
for a nucleus to develop a halo is that the outmost nucle- 
ons move in shells which extend far in space. That is, only 
a weak barrier keep the system within the nuclear vol- 
ume. These shells may be resonances, antibound states 
(also called virtual states), or even low-spin bound states 
which lie very close to the continuum threshold. These 
conditions are fulfilled by the nucleus 11 Li and also heav- 
ier Li isotopes. There are a number of experiments which 
have been performed in these very unstable isotopes in or- 
der to get information about the structure of halos [I[ . In 
particular we will concentrate our attention to Refs. 0- 
[7[ where the spectrum of 12 Li was measured. Our aim 
is to analyze these experimental data by using a suitable 
formalism to treat unstable nuclei. This formalism is an 
extension of the shell model to the complex energy plane 
and is therefore called complex shell model 8j, although 
the name Gamow shell model is also used 9|. In addi- 
tion, the correlations induced by the pairing force acting 
upon particles moving in decaying single-particle states 
will be taken into account by using the multistep shell 
model (MSM) 0. 

The formalism is presented in Section [XT] Applications 
are in Section Mil and a summary and conclusions are in 
Section HVl 



II. THE FORMALISM 

The study of unstable nuclei is a very difficult un- 
dertaking since, in principle, time dependent formalisms 



should be used to describe the motion of a decaying nu- 
cleus. However, the system may be considered stationary 
if it lives a long time. In this case the time dependence 
can be circumvented. In fact, often unstable nuclei live 
a very long time and therefore they may be considered 
bound as, e.g., in alpha decaying states of many heavy 
isotopes, like 208 Bi or 180 Ta(9-), with T 1/2 > 10 15 y. On 
the other hand, experimental facilities allow one nowa- 
days to measure systems living a very short time. To 
describe these short time processes one has to consider 
the decaying character of the system. This is shown, e.g ., 
by the failure of the shell model calculation of Ref. [Lj| , 
performed by using a standard bound representation, to 
explain even the ground state of 12 Li [5[- 

Of the various theories that have been conceived to 
analyze unbound systems, we will apply an extension of 
the shell model to the complex energy plane Q . The ba- 
sic assumption of this theory is that resonances can be 
described in terms of states lying in the complex energy 
plane. The real parts of the corresponding energies are 
the positions of the resonances while the imaginary parts 
are minus twice the corresponding widths, as it was pro- 
posed by Gamow at the beginning of quantum mechan- 
ics [Hj] • These complex states correspond to solutions of 
the Schrodinger equation with outgoing boundary con- 
ditions. We will not present here the formalism in de- 
tail, since this was done many times before, e.g., in Refs. 
[ill, [13]. Rather, we will give the main points necessary 
for the presentation of the applications. 



A. The Berggren representation 

In this Subsection we will very briefly describe the rep- 
resentation to be used here. 

The eigenstates of a central potential obtained as out- 
going solutions of the Schrodinger equation can be used 
to express the Dirac 5-function as |15| 



(J(r-r') = y^ y w„(r)w n (r')+ / dEu{r,E) 



u(r',E), (1) 
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where the sum runs over all the bound and antibound 
states plus the complex states (resonances) which lie be- 
tween the real energy axis and the integration contour 
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L + . The wave function of a state n in these discrete set 
is w n (r) and u(r,E) is the scattering function at energy 
E. The antibound states are virtual states with nega- 
tive scattering length. They are fundamental to describe 
nuclei in the Li region (l6j . 

Discretizing the integral of Eq. (JlJ one obtains the set 
of orthonormal vectors forming the Berggren repre- 
sentation [17] . Since this discretization provides an ap- 
proximate value of the integral, the Berggren vectors ful- 
fill the relation / w . \ipj)(ipj |, where all states, that is 
bound, antibound, resonances and discretized scattering 
states, are included. The corresponding single-particle 
wave functions are 

W<Pi)=RTHid t (r)(xi/Mr)) jimi , (2) 

where x is the spin wave function and 

Rmhjiir) = <t>mhji(r)/r (3) 

is the radial wave function fulfilling the Berggren metric, 
according to which the scalar product between two func- 
tions consists of one function times the other (for details 
see Ref. [13), i- e -, 



(4) 



Using the Berggren representation one readily gets the 
two-particle shell-model equations in the complex energy 
plane (CXSM) [3, i.e., 



(W(a 2 )-e l -e J )X(ij; a 2 ) = ^{kl; a 2 \V\ij; a 2 )X(kl; a 2 ), 

k<l 

(5) 

where V is the residual interaction. The tilde in the in- 
teraction matrix element denotes mirror states so that in 
the corresponding radial integral there is not any com- 
plex conjugate, as required by the Berggren metric. The 
two-particle states are labeled by a 2 and Latin letters 
label single-particle states. W(a 2 ) is the correlated two- 
particle energy and e$ is single-particle energy. The two- 
particle wave function is given by 

\a 2 )=P+(a 2 )\0), (6) 

where the two-particle creation operator is given by 



P + (a 2 ) 



X(ij;a 2 )- 



(7) 



and A Q2 is the angular momentum of the two-particle 
state. 

We will use a separable interaction as in Ref. fl8j], 
which describes well the states of 11 Li. The energies are 
thus obtained by solving the corresponding dispersion re- 
lation. The two-particle wave function amplitudes are 
given by 



X(ij;a 2 ) = N ai 



where f(ij,ct 2 ) is the single particle matrix element of 
the field defining the separable interaction and N a2 is 
the normalization constant determined by the condition 

The spectrum of 11 Li was already evaluated within the 
CXSM including antibound states [Hj]. Here we will 
repeat that calculation in order to determine the two- 
particle states to be used in the calculation of the three- 
particle system, i.e., 12 Li. For this we will use the Multi- 
step Shell Model Method, which we will briefly describe 
below. 



B. The Multistep Shell Model Method 

As its name indicates, the Multistep Shell Model 
Method (MSM) solves the shell model equations in sev- 
eral steps. In the first step the single-particle representa- 
tion is chosen. In the second step the energies and wave 
functions of the two-particle system are evaluated by us- 
ing a given two-particle interaction. The three-particle 
states are evaluated in terms of a basis consisting of the 
tensorial product of the one- and two-particle states pre- 
viously obtained. In this and subsequent steps the inter- 
action does not appear explicitly in the formalism. In- 
stead, it is the wave functions and energies of the compo- 
nents of the MSM basis that replace the interaction. The 
MSM basis is overcomplete and non-orthogonal. To cor- 
rect this one needs to evaluate the overlap matrix among 
the basis states also. A general description of the formal- 
ism is in Ref. fl9l ] . The particular system that is of our 
interest here, i.e., the three-particle case, can be found 
in Ref. [l(|, where the MSM was applied to study the 
three-neutron hole states in the nucleus 205 Pb. 

Using the Berggren single-particle representation de- 
scribed above, we will evaluate the complex energies and 
wave functions of 12 Li using the MSM basis states con- 
sisting of the Berggren one-particular states, which are 
states in 10 Li, times the two-particle states corresponding 
to 11 Li. Below we refer to this formalism as CXMSM. 

The three-particle energies W(a 3 ) are given by [l(3] 

{W(a 3 ) - Ei - W{a 2 ))(a 3 \(c+P+{a 2 )) a3 \0) 
= E I E - e « - £ k )A(ia 2 ,jp 2 ; k) 1 



x (a 3 \(cfP+(l3 2 )) a3 \0), (9) 



where 



A(ia 2 ,j/3 2 ;k) =a 2 fo{ * * f 2 \ Y(kj; a 2 )Y(ki;j3 2 ), 



J a 3 a 2 



(10) 



(8) 



Y(ij;a 2 ) = {l + 5{i,j)) 1/2 X[iy,a 2 ), (11) 
and the rest of the notation is standard. 
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The matrix defined in Eq. © is not hermitian and the 
dimension may be larger than the corresponding shell- 
model dimension. This is due to the violations of the 
Pauli principle as well as overcounting of states in the 
CXMSM basis. Therefore the direct diagonalization of 
Eq. (J9]) is not convenient. One needs to calculate the 
overlap matrix in order to transform the CXMSM basis 
into an ortho normal set. In this three-particle case the 
overlap matrix is 

(0|(c+P + M)i 3 (c+P+(/3 2 )) Q3 |0) 

= SijS a2 2 + y ^2,A(ia.2,jp2)k). (12) 

k 

Using this matrix (fT2"|) one can transform the ma- 
trix determined by Eq. Q into a hermitian matrix T 
which has the right dimension. The diagonalization of 
T provides the three-particle energies. The correspond- 
ing wave function amplitudes can be readily evaluated to 
obtain 

\a 3 }=P + (a 3 )\0), (13) 
P + (a 3 ) = Y / X(^2;a 3 )(c+P+(a 2 )) a3 , (14) 

JCK2 

where P + (a 3 ) is the three-particle creation operator. 

It has to be pointed out that in cases where the basis 
is overcomplete the amplitudes X are not well defined. 
But this is no hinder to evaluate the physical quantities. 
For details see Ref. [13 ■ 

The advantage of the MSM in stable nuclei is that one 
can study the influence of collective vibrations upon nu- 
clear spectra within the framework of the shell model. 
Thus, in Ref. [l(| the multiple structure of particle- 
vibration coupled states in odd Pb isotopes was analyzed. 

But the most important feature for our purpose is that 
the CXMSM allows one to choose in the basis states a 
limited number of excitations. This is because in the 
continuum the vast majority of basis states consists of 
scattering functions. These do not affect greatly phys- 
ically meaningful two-particle states. That is, the ma- 
jority of the two-particle states provided by the CXSM 
are complex states which form a part of the continuum 
background. Only a few of those calculated states cor- 
respond to physically meaningful resonances, i.e., reso- 
nances which can be observed. Below we call a "res- 
onance" only to a complex state which is meaningful. 
These resonances are mainly built upon single-particle 
states which are either bound or narrow resonances. Yet, 
one cannot ignore the continuum when evaluating the 
resonances. The continuum configurations in the reso- 
nance wave function are small but many, and they af- 
fect the two-particle resonance significantly [l3[. That 
is, the important continuum configurations to induce res- 
onances are contained in the corresponding resonance 
wave functions. Therefore, the great advantage of the 
CXMSM is that one can include in the basis only two- 
particle resonances, while neglecting the background con- 
tinuum states, which form the vast majority of complex 



two-particle states. The question is which are the two- 
particle states that are indeed resonances. This, and also 
the question of how to evaluate and recognize the phys- 
ically meaningful three-particle states, are addressed in 
the next Section. 



III. APPLICATIONS 

In this Section we will apply the CXMSM formalism 
described above to study the nucleus 12 Li. As in Ref. 
[lH, we will take the core to be the nucleus 9 Li. This 
is justified since the three protons included in this core 
can be considered frozen and, therefore, merely specta- 
tors 0. 

To evaluate the valence shells we will proceed as in 

Refs. mm and choose as central field a Woods-Saxon 

potential with different depths for even and odd orbital 
angular momenta I. The corresponding parameters are 
(in parenthesis for odd Z- values) a— 0.670 fm, r = 1.27 
fm, Vn = 50.0 (36.9) MeV and V so =16.5 (12.6) MeV. As in 
Ref. [23|, we thus found the single-particle bound states 
Osi/2 at -23.280 MeV and 0p 3/2 at -2.589 MeV forming 
the 9 Li core. The valence shells are the low lying reso- 
nances 0pi/2 at (0.195,-0.047) MeV and 0d 5/2 at (2.731, 
-0.545) MeV and the shell 0d 3/2 at (6.458,-5.003) MeV. 
This cannot be considered as a resonance, since it is so 
wide that rather it is a part of the background. Besides, 
the state lsi/ 2 appears as an antibound state at -0.050 
MeV. We thus reproduce the experimental single-particle 
energies as given in Ref. [24j . We also found other states 
at higher energies, but they do not affect our calculation 
because they are very high and also very wide to be con- 
sidered as meaningful resonances. We thus include in our 
Berggren representation only the antibound state lsi/2 
and the resonances 0pi/ 2 and 0d 5 / 2 . 

The energy of the resonance 0p\/ 2 has been contested 
and instead the value (0.563,-0.252) MeV was proposed 
[IJ. Since this is an important quantity in the determi- 
nation of the two- and three-particle spectra, we will use 
both of them in our calculations. We obtained this 0p\/ 2 
energy (i.e., (0.563,-0.252) MeV) by choosing V =34.755 
MeV for I odd while keeping all the other parameters as 
before. 

To define the Berggren single-particle representation 
we still have to choose the integration contour L + (sec 
Eq. ©). 

To include in the representation the antibound lsi/2 
state as well as the Gamow resonances 0pi/ 2 and 0d 5 / 2 we 
will use two different contours. The number of points on 
each contour defines the energies of the scattering func- 
tions in the Berggren representation, i.e., the number of 
basis states corresponding to the continuum background. 
This number is not uniformity distributed, since in seg- 
ments of the contour which are close to the antibound 
state or to a resonance the scattering functions increase 
strongly. We therefore chose the density of points to be 
larger in those segments. 
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We include the antibound state by using the contour 
in Fig. [T] 



With the single-particle representation thus dehned we 
proceed to evaluate the two-particle states. 
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FIG. 1. Contour used to include the antibound state (see, 
also, Ref. [13]). The points Bi denote bound state energies 
while A denotes the antiboud state. The points Vt corre- 
spond to the vertices defining the contour. They have the 
values Fi=(-0.05,0.05) MeV, F 2 =(-0.1,0) MeV, F 3 =(0,-0.4) 
MeV, F 4 =(0.5,-0.4) MeV, F 5 =(8,-0.4) MeV, V 6 =(8,0) MeV 
and W=(10,0) MeV 

The number of points in each segment are given in 
Table H 




Re(E) 



FIG. 2. Contour used to include the Gamow resonances 
represented by the points Gi. The vertices are Vi = (0,-1) 
MeV, V2=(l-1) MeV, V 3 =(8,-l) MeV, Vi=(8,0) MeV and 
V 5 =(10,Q) MeV. 

For the Gamow resonances the contour in Fig. [5] is 
used with the number of Gaussian points as in Table |TTJ 

We have adopted these points after verifying that the 
results converged to their final values. A discussion about 
the choice of these contours and also on the physical 
meaning of the antibound state can be found in Ref. [l8[ . 



A. Two- particle states: the nucleus 11 Li. 

The only state which is measured in 11 Li is its bound 
ground state, which was found to lie at an energy of - 
0.295 MeV [25]. However, this state is more bound than 
that, as it was determined in more recent experiments 
[26M28I ] . We will adopt the most precise of these val- 
ues, i.e., -0.369 MeV [27|. The corresponding angular 
momentum is 3/2~. This spin arises from the odd pro- 
ton lying deep in the spectrum coupled to two neutrons. 
As has been already mentioned, the proton is considered 
to be a spectator. The dynamics of the system is thus 
determined by the pairing force acting upon the two neu- 
trons coupled to a state + .which behaves as a normal 
even-even ground state [2l[ . Besides the energy, this 
state has been measured to have an angular momentum 
contain of about 60 % of s- waves and 40 % of p-waves, al- 
though small components of other angular momenta are 
not excluded (25[. 

We will perform the calculation of the two-particle 
states by using the separable interaction discussed in Sec- 
tion [TTJ The strength G\ 2 , corresponding to the states 
with angular momentum A2 and parity (— l) Aa , will be 
determined by fitting the experimental energy of the low- 
est of these states, as usual. It is worthwhile to point 
out that G\ 2 defines the Hamiltonian and, therefore, is 
a real quantity. The two-particle energies are found by 
solving the corresponding dispersion relation while the 
two-particle wave function components are as in Eq. ([H]). 

The angular momentum contains of the ground state 
wave function are shown in Table IIII1 One sees that 
for the case in which the single-particle state 0pi/ 2 is 
assumed to lie at (0.195,-0.047) MeV the two-particle 
wave function consists of 46.8% s-states and 49.1% p- 
states which are reasonable values. For the Opi/2 energy 
at (0.563,-0.252) MeV the angular momentum contain is 
72.6 % s-states and 20.9 % p-states, which is also accept- 
able, specially considering that it provides the correct 
order of the relative magnitudes. Both cases are in rea- 
sonable agreement with experiment. 

The wave function components corresponding to this 
state are strongly dependent upon the contour that one 
uses. However, measurable quantities, like the energies 
and transition probabilities, do not. This is because the 
physical quantities are defined on the real energy axis 
and, therefore, they remain the same when changing con- 
tour. But complex states which are part of the continuum 
background do not have any counterpart on the real en- 
ergy axis and the physical quantities for these states ac- 
quire different values for different contours [lq |. We will 
use this property to determine whether a complex state 
is a meaningful resonance. This is important, since the 
ground state is the only one for which experimental data 
exists. There might be other meaningful states that have 
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TABLE I. Number of Gaussian points in the different segments of the contour of Fig. [TJ 



Segment 


[(0,0) -Vi] 


[Vi - v 2 ] 


[Vi - V 3 ] 


[V 3 - Vi] 


[Vi - V B ] 


[V 5 - V 6 ] 


We - Vr] 


Number 


30 


30 


30 


30 


30 


16 


6 



TABLE II. Number of Gaussian points in the different seg- 



ments of the contour of Fig. [2] 


Segment [(0, 0) - Vi] [Vi - V 2 ] [V 2 - Va] [V 3 - Vi] 
Number 30 30 30 8 


[Vi - Vs] 
4 


TABLE III. Angular momentum contain of 11 Li(f 
sponding to the energies e Pl , 2 discussed in the text. 


;s) corre- 


component(%) 


e Pl / a (MeV) s-waves p-waves 


d-waves 


(0.195,-0.047) 46.8 49.1 
(0.563,-0.252) 72.6 20.9 


4.2 
6.4 



not been found yet. This implies that we have to evaluate 
all possible two-particle states which are spanned by our 
single-particle representation. To decide whether a state 
thus calculated is a meaningful resonance we will pro- 
ceed as in Refs. 0, H^] and analyze the singlet (5 = 0) 
component of the two-particle wave function. The cor- 
responding expression for this component was given in 
Eq. (10) of Ref. [29j |. but we will show it here again for 
clarity of presentation. For the state a with spin and 
spin-projection (JM) that component is, with standard 
notation, 

tfajM^) = [xi/2(l)Xi/2(2)]2$Z X(ab,aJM)j a j b 

a<b 

x [C(ab,nr 2 ) - (-)^+^- J C(&a,r 1 r 2 )] (15) 

where 

C(a6,fir 2 ) = Mn)Mr2)H h+1/2 - la+J 

x \la 3 a l /A [Yla {h)Y lb {f 2 )] 1 :^ (16) 

[Jbh J J 

and <p a {r) is the radial wave function corresponding to 
the single-particle state a (Eq. j3j). 

If the two-particle state (aJM) is a meaningful reso- 
nance then the wave function above should be localized 
within a region extending not too far outside the nuclear 
surface, and its imaginary part should not be too large 
18]. To study these features it is not necessary to go 
to all six dimensions corresponding to the coordinates r\ 
and r*2. In fact it is enough to consider the coordinate 
r given by f i = r- 2 = f — (0, 0, r) which corresponds to 
the two particles located at the same point and in the 



z-direction. For details see 29] . We will call this one- 
dimensional function ^ a jAi(r)- 

The evaluation of the + states is a relatively easy task, 
since in this case we have determined the strength G + 
by fitting the experimental energy of 11 Li(gs). With this 
value of the strength we calculated all the + states and 
found that there is no any meaningful resonance. The 
only physically meaningful + state is the bound ground 
state. Besides this, we found that there is a meaningful 
resonance, which is the state 2+ at an energy (2.300,- 
0.372) MeV. In Fig. [5] we show the corresponding radial 
wave function ^f 2 + (r) . One sees that it is rather local- 
ized and that its imaginary part is relatively small as 
compared to the corresponding real part. This is a state 
which perhaps is at the limit of what can be considered 
a meaningful resonance. Yet, it is not too wide and it 
has an effect on the physical three-particle states, as will 
be seen below. It is worthwhile to point out that the 
width of this state (744 keV) is the escape width. At 
high energies, where the giant resonances lie, most of the 
width consists of the spreading width, i.e., of mixing with 
particle-hole configurations 30]. However, at the low en- 
ergies of the states that we study this mixing would not 
be relevant. 

We have thus found that the only two-particle states 
on which physically meaningful three-particle states can 
be built are 11 Li(gs) and n Li(2^). One can understand 
why there are so few meaningful two-particle states in 
this nucleus by looking at the radial wave functions of 
the single-particle states that form the representation. 
For the antibound state one sees in Fig. [3] that it ex- 
tends in an increasing rate far out from the nuclear sur- 
face, as expected in this halo nucleus (the standard value 
of the radius is here 1.2 x ll 1 / 3 =2.7 fm). The radial wave 
function corresponding to the Gamow resonance Opi / 2 is 
shown in Fig. 2] The state n Li(2^) is determined by 
the antibound state and the resonance, which has 

a large and increasing imaginary part at relative short 
distances, as shown in Fig. [5] These single-particle wave 
functions have large imaginary parts and a divergent na- 
ture even at rather short distances from the nucleus. In 
other words, they correspond to states that live a time 
which is too short to produce meaningful two-particle 
resonances. 

An important point for the analysis of the three- 
particle states to be performed below, is that the scatter- 
ing wave functions in the segments [(0, 0) — Vi], [Vi — V2] 
and [V 2 — V3] are similar in magnitude to the wave func- 
tion of the antibound state. This is because the segments 
are very close to the antibound state 18]. This is a fea- 
ture that cannot be avoided, and is due to the attractive 



6 



character of the pairing force. That is, the lowest single- 
particle configuration in the Berggren basis is V^ 2 , with 
energy — 2e, where V2 = (— e,0). This configuration has 
to lie above the energy of the two-particle correlated state, 
i. e. it has to be e > w( 11 Li(gs))/2. 
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FIG. 5. As Fig. [3] but for the Gamow resonance Qd 5 / 2 at an 
energy of (2.731,-0.545) MeV. The dashed line is the imagi- 
nary part of the wave function. 



FIG. 3. Radial function <j>(r) corresponding to the single- 
particle neutron anti bound state Osj/2 at an energy of -0.050 
MeV. 




10 12 14 



r[fm] 



FIG. 4. As Fig. [3] but for the Gamow resonance Opi/2 at an 
energy of (0.563,-0.252) MeV. The dashed line is the imagi- 
nary part of the wave function. 

With the states 0^ and 2^ thus calculated we pro- 
ceeded to the calculation of the three-particle system 
within the CXMSM. 



B. Three-particle states: the nucleus 1 Li. 

With the single-particle states and the two-particle 
states 11 Li(gs) and n Li(2+) discussed above, we formed 




FIG. 6. Radial function \l/(r) corresponding to the two- 
particle state U Li(2+) at an energy of (2. 300,-0. 372)MeV. The 
dashed line is the imaginary part of the wave function. 



all the possible three-particle basis states. We found that 
the only physically relevant states are those which are 
mainly determined by the bound state 11 Li(0 ] t "). The cor- 
responding spins and parities are 1/2+, 1/2 - and 5/2 + . 
States like 3/2 + , which arises from the CXMSM config- 
uration I lsi/2 <8> 2^; 3/2 + >, is not a meaningful state. 

Due to the large number of scattering states included in 
the single-particle representation the dimension of three- 
particle basis is also large. The scattering states are 
needed in order to describe these unstable states. 

In the calculations we took into account all the pos- 
sibilities described above regarding the energies of the 
single-particle state Opi/2 as well as the binding energy of 
the state n Li(gs). The corresponding results are shown 
in Tables HVl and fVl Comparing these two Tables one sees 
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TABLE IV. Calculated three-particle states in 12 Li (in MeV) 
corresponding to the two energies e Pl , 2 of Table IIIII and the 
two-particle energies taken to be u Q + — — 0.295MeV and 

u) 2 + = (2.300, -0.372)MeV. 



e Pl/2 


1/2+ 


1/2- 


5/2+ 


(0.195,-0.047) 
(0.563,-0.252) 


(-0.386,-0.006) 
(-0.382,-0.006) 


(0.821,-0.189) 
(1.114,-0.403) 


(1.348,-0.276) 
(1.169,-0.242) 


TABLE V. 


As Table |IV] but for w 0+ = - 


0.369MeV. 




1/2+ 


1/2- 


5/2+ 


(0.195,-0.047) 
(0.563,-0.252) 


(-0.466,-0.011) 
(-0.466,-0.011) 


(0.753,-0.206) 
(1.116,-0.411) 


(1.315,-0.276) 
(1.148,-0.243) 



that the state 1/2+ depends very slightly on the energy 
e Pi/2 but rather strongly on uj a . 2 . Instead, this tendency 
is opposite for the states l/2~ and 5/2+. 

The most important feature in these Tables is that the 
lowest state is 1/2+ and that it has a real and negative 
energy. It is an antibound state, as it is the lsi/2 state 
itself. A manifestation of this is that the radial wave 
function diverges at large distances. 

Accepting the latest reported values for the energies 
of the states Opi/2 and 11 Li(gs), i.e., those in the second 
line of Table [V] our calculation predicts, besides the an- 
tibound 1/2+ state, a state 1/2" at (1.116,-0.411) MeV 
and a state 5/2+ at (1.148,-0.243) MeV. This assignment 
agrees well with what is given in Ref. Q for the ground 
state of 12 Li, which was found to be an antibound (or 
virtual) state. In the same fashion, in Ref. a state 
was found at 1.5 MeV which probably has spin and par- 
ity 5/2+. In these experiments the angular momentum 
contain of the states were measured and therefore it is 
proper to compare the experimental quantities with our 
calculations, where only neutrons are considered. 

In Ref. [3] it was also found that 12 Li(gs) is an anti- 
bound state but, in addition, two other low-lying states 
were observed at 0.250 MeV and 0.555 MeV by using two- 
proton removal reactions. In this case the 0j?3/ 2 proton in 
the core may interfere with the neutron excitations eval- 
uated above. In particular, the antibound 1/2+ ground 
state would provide, through the proton excitation, a 
state 1~ and a 2". This is the situation encountered in 
the shell model calculation presented in Ref. |7j. As we 
have pointed out above, from the CXMSM viewpoint the 
very unstable states determining the spectrum in this nu- 
cleus, with wave functions which are both diverging and 
complex, can hardly be described by harmonic oscillator 
representations. To investigate the relation between the 
results provided by both representations we repeated the 
shell-model calculation of Ref. Q- We thus took 4 He as 
the core and used the WBT effective interaction [3l|, [HJ ■ 
The resulting Hamiltonian matrix was diagonalized by 



using the code described in 33] . The corresponding cal- 
culated energies, which agree with those presented in Q, 
are shown in Fig. [7] 

One sees in this Figure that the full calculation pre- 
dicts all excited states to lie well above the correspond- 
ing experimental values. It is worthwhile to point out 
that the calculated states exhibit rather pure shell model 
configurations. For instance the states 2^ (ground 
state) and 1^ are mainly composed of the configuration 
k [OP3/2] v [(0pi/ 2 ) 2 lsi /2 ]). This does not fully agree 
with our CXMSM calculation, since in our case this wave 
function is mainly of the form |lsi/2 ® n Li(gs)). This 
differs from the shell model case in two ways. First, 
the state 11 Li(gs) contains nearly as much of lsi/ 2 as of 
0pi/2- Second the continuum states contribute much in 
the building up of the antibound 12 Li(gs) wave function, 
as discussed above. In our representation it is straight- 
forward to discern the antibound character of this state, 
which is not the case when using harmonic oscillator 
bases. 

The shell model splittings between states originated 
from the same configurations in Fig. [7] are due to 
the neutron-proton interaction, which in some cases 
can be large. For instance, the matrix elements 
(7r[0p 3/2 ]i/[0£i /2 ]; J|V|7r[0p3/ 2 ]^[0pi/ 2 ]; J) are strong and 
attractive [34|. 

Given all the uncertainties related to these calculations 
in the continuum, and the scarce amount of experimental 
data, we will not attempt to evaluate the 12 Li states by 
adding a new uncertainty which would be the inclusion 
of a proton-neutron interaction. 

Considering the limitations that one expects from a 
calculation of very unstable states by using harmonic 
oscillator basis, the rather good agreement between the 
shell-model and the CXMSM presented above is surpris- 
ing. In contrast to the shell-model, the CXMSM pro- 
vides not only the energies but also the widths of the 
resonances. It is therefore very important to encourage 
experimental groups to try to obtain these quantities in 
order to probe the formalisms. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have studied excitations occurring in 
the continuum part of the nuclear spectrum which are 
at the limit of what can be observed within present ex- 
perimental facilities. These states are very unstable but 
yet live a time long enough to be amenable to be treated 
within stationary formalisms. We have thus adopted the 
CXSM (shell model in the complex energy plane [l3| ) for 
this purpose. In addition we performed the shell model 
calculation by using the multistep shell model. In this 
method of solving the shell model equations one proceeds 
in several steps. In each step one constructs building 
blocks to be used in future steps (l9| . We applied this 
formalism to analyze 12 Li as determined by the neutron 
degrees of freedom, i.e., the three protons in the core 
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FIG. 7. Experimental level scheme in 12 Li. The three lowest levels are from [7J, while the one at 1.5 MeV is from |fj]. In the 
second column are the three-neutron CXMSM results. In the columns A-C are the shell-model calculations corresponding to 
different truncation schemes: A) within — lfku excitations, B) within — 3fko excitations and C) full psd space. Dashed lines 
indicate the widths of the resonances. 



were considered to be frozen. In this case the excitations 
correspond to the motion of three particles, partitioned 
as the one-particle times the two-particle systems. This 
formalism was applied before, e.g., to study multiplets in 
the lead region flOf . 

By using single-particle energies (i.e., states in 10 Li) as 
provided by experimental data when available or as pro- 
vided by our calculation, we found that the only phys- 
ically meaningful two-particle states are 11 Li(gs), which 
is a bound state, and 11 Li(2^), which is a resonance. As 
a result there are only three physically meaningful states 
in 12 Li which, besides the antibound ground state, it is 
predicted that there is a resonance l/2~ lying at about 
1 MeV and about 800 keV wide and another resonance 
which is 5/2 + lying at about 1.1 MeV and 500 keV wide. 
That the ground state is an antibound (or virtual) state 
was confirmed by a number of experiments J5j-l7| and the 
state 5/2+ has probably been observed in @. However, 



in [7J two additional states, lying at rather low energies, 
have been observed which do not seem to correspond to 
the calculated levels. It has to be mentioned that neither 
a shell model calculation, performed within an harmonic 
oscillator basis, provides satisfactory results in this case. 
Yet, we found that this shell model calculation works 
better than one would assume given the unstable char- 
acter of the states involved. We therefore conclude that 
in order to probe the formalisms that have been intro- 
duced to describe these very unstable systems additional 
experimental data, specially regarding the widths of the 
resonances, would be required. 
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